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Charm and bottom quark masses from QCD moment sum rules 
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In this work the charm and bottom quark masses are determined from QCD moment sum rules for the charmonium 
and upsilon systems. In our analysis we include both the results from non-relativistic QCD and perturbation theory 
at next-next-to-leading order. For the pole masses we obtain M c ~ 1.75 ± 0.15 GeV and Mb = 4.98 ± 0.125 GeV. 
Using the potential-subtracted mass in intermediate steps of the calculation the MS-masses are determined to 
m c (m c ) = 1.19 ± 0.11 GeV and m b (m b ) = 4.24 ± 0.10 GeV. 



1. Introduction 

An important task within modern particle phe- 
nomenology consists in the determination of the 
quark masses, being fundamental parameters of 
the Standard Model. In the past, QCD mo- 
ment sum rule analyses have been successfully ap- 
plied for extracting the charm and bottom quark 
masses from experimental data on the charmo- 
nium and bottomium systems respectively @,^,^| . 
The basic quantity in these investigations is the 
vacuum polarisation function II(g 2 ): 



n 



»Al 2 ) = * J d 4 xe^(T{j M (x)jt(0)}) 



(i) 



where the relevant vector current is represented 
either by the charm j^(x) = {c^^c){x) or the 
bottom current j b Ax) — (bj^b) (x). Via the 
optical theorem, the experimental cross section 
a(e + e~ — > cc, bb) is related to the imaginary part 
of 11(a): 



R(s) = 



1 er(e + e 



cc, bb) 



= 127rimII(s) . 



Usually moments of the vacuum polarisation are 
defined by taking derivatives of the correlator at 
s = 0. However, in this work we allow for an 
arbitrary evaluation point s = — 4m 2 £ to define 
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the dimensionless moments 



12tt 2 



4m 



ds 



IL(s) 



s — — 4m 2 £ 



- 2 l±- v{1 - v2)n " R{v) (2) 



where v = i/T- 4to 2 /s is the velocity of the 
heavy quark. The parameter £ encodes much in- 
formation about the system. By taking £ larger 
the evaluation point moves further away from the 
threshold region. Consequently, the theoretical 
expansions show a better convergence, but at the 
same time the sensitivity on the mass is reduced. 
The moments can either be calculated theoret- 
ically, including Coulomb resummation, pertur- 
bation theory and nonperturbative contributions, 
or be obtained from experiment. In this way one 
can relate the heavy quark masses to the hadronic 
properties of the quark- antiquark systems. 

A characteristic feature of these heavy-heavy 
bound state systems is the Coulomb-like form 
of the potential. Developing the quantum me- 
chanical sum rules for Coulomb systems, in jj] 
it has been shown that the application of fixed- 
order perturbation theory in such systems results 
in unstable predictions for the masses. To ob- 
tain a stable sum rule it is necessary to incorpo- 
rate the threshold behaviour which can be calcu- 
lated in the framework of non-relativistic QCD 
(NRQCD). 
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A natural choice for the mass in eq. (g) is the 
pole mass M. First we will employ the pole mass 
scheme to extract the pole masses. However, as 
the pole masses suffer from renormalon ambigui- 
ties H , we shall then use the potential-subtracted 
(PS-) mass mps f(| to extract the MS-masses. 

2. Coulomb resummation 

The theory of NRQCD provides a consistent 
framework to treat the problem of heavy quark- 
antiquark production close to threshold. The 
contributions can be described by a nonrelativis- 
tic Schrodinger equation and systematically cal- 
culated in time-independent perturbation theory 
(TIPT) . The correlator is expressed in terms of a 
Green's function G(k) = G(0,0,fc) @,|,§: 



H( S ) 



2M 2 



Ak 2 

C h (a s )G(k) + ^Gc(fc) ) , (3) 



where k = -J M 2 — s/4 and M represents the pole 
mass. The constant Ch{a s ) is a perturbative co- 
efficient needed for the matching between the full 
and the nonrelativistic theory and naturally de- 
pends on the hard scale. 

To calculate the moments from the Green's 
function we will directly perform the derivatives 
at s = — 4M 2 £ according to eq. (||). Since the 
Green's function is known analytically Q as a 
function of k = k{s), this can be done numeri- 
cally. In this way we take advantage of the fact 
that the perturbative expansion parameter de- 
pends on the evaluation point. In addition, we 
can extract the spectral density above threshold 
by taking the imaginary part of the Green's func- 

tion gg. 

The moments depend on the three scales Hsoft, 
l-ifac an d fJ-hard, the soft, factorisation and hard 
scale respectively. The residual dependence on 
these scales will turn out to give the dominant er- 
ror in the determination of the masses. The large 
corrections are partly due to the definition of the 
pole mass. These contributions can be reduced 
by using an intermediate mass definition. In this 
analysis we will use the potential-subtracted (PS) 
mass H where the potential below a separation 
scale fi sep is subtracted. This mass definition 
leads to an improved scale dependence and a more 



precise determination of the MS-masses. 



3. Perturbative expansion 

The perturbative spectral function R Pert (s) 
can be expanded in powers of the strong coupling 
constant, 

R pt ( S ) = i?(°)(s) + 2- i?(D( s ) + ^ i? (2) (s) + • • ■ 



From this expression the corresponding moments 
M. n can be calculated via the integral of eq. (||). 
The first two terms are known analytically and 
can for example be found in ref. |11|. R^ 2 \s) 
is still not fully known analytically. We employ a 
method based on Pade-approximants to construct 
the spectral density in the full energy range fl2] , 
It uses available information around q 2 = 0, 
at threshold and in the high energy region. It 
has the advantage that it gives a good description 
until relatively close to threshold, the region on 
which the quark masses are most sensitive. 



4. Condensate contributions 

The non-perturbative effects on the vacuum 
correlator are parametrised by the condensates. 
The leading correction is the gluon condensate 
contribution which is known up to next-to-leading 
order jl4|. Furthermore, in ]l5|,^6[ the dimension 
6 and 8 condensates have been calculated. From 
the numerical analysis it turns out that the ab- 
solute contribution of the condensates to the mo- 
ments is negligible for the bottomium and of little 
influence for the charmonium. The relative sup- 
pression of the gluon condensate to former char- 
monium analyses is due to three reasons: First, 
the absolute value of the theoretical moments in- 
creases from the Coulomb resummation. Then we 
evaluate the moments at larger £ and smaller n 
where the nonperturbative contributions are rel- 
atively small. Finally, since we obtain a larger 
pole mass than former analyses, the condensates, 
starting with a power of 1/M 4 , are suppressed 
further. 
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5. Phenomenological spectral function 



4.. 



Experimentally, the six lowest lying ip— and 
T— resonances have been observed. Furthermore, 
recent measurements of BES [[l7| in the charmo- 
nium region have improved the cross section be- 
tween 3.7 and 4.8 GeV. Since the widths of the 
poles are very small compared to the masses, the 
narrow-width approximation provides an excel- 
lent description of these states. To model the 
contributions above the 6th resonance in the bot- 
tomium system we use the assumption of quark- 
hadron-duality and integrate the theoretical spec- 
tral density above ^/Iq = 11.0 ± 0.2 GeV. In 
the charmonium system we include the two low- 
est resonances, the BES data and the theoretical 
spectral density above 4.8 GeV. 

6. Numerical analysis 

The theoretical part of the correlator contains 
the poles of the Green's function, the spectral 
density above threshold and the condensates. For 
high velocities the spectral density is well de- 
scribed by the perturbative expansion whereas 
the resummed spectral density gives a good de- 
scription for low values of v. Therefore we con- 
struct a theoretical spectral density in the full 
energy range which includes the essential infor- 
mation in both regions of v. For a more detailed 
discussion the reader is referred to Now we 
discuss the most important points in the numer- 
ical analysis of the bottomium and charmonium 
systems respectively. 

6.1. Bottom mass 

First we discuss the values of £ and n. Since 
the bottom quark is relatively heavy, even for 
£ = the nonrelativistic and perturbative ex- 
pansions converge reasonably well. Nevertheless, 
the contributions from the poles of the Green's 
function still dominate the theoretical part. To 
reduce their influence and to spread the theoret- 
ical contributions more equally among the poles, 
the resummed spectral density and the perturba- 
tive spectral density we must choose a higher £. 
However, for £ > 1 the moments loose sensitivity 
on the mass and the error from the input param- 
eters increases. Therefore we use a central value 
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Figure 1. Thick solid line: central PS-mass; thin 
solid lines: mps^ for fX S oft — 2.0 and 3.5 GeV; 
dashed lines: mps,& for /i/ ac = 2.0 and 5.0 GeV; 
dotted lines: mps j0 for Hhard — 2.5 and 10.0 GeV. 

of £ = 0.5 and vary £ between < £ < 1. Since 
the relevant scale for the evaluation point is the 
lowest bound state energy, values of £ = 0, 0.5 or 
1 already correspond to well separated evaluation 
points. For n we use a range of 5 < n < 10 where 
the theoretical expansion and the phenomenolog- 
ical uncertainty are under control. 

As central values for our scales we have se- 
lected fisoft = 2.5 GeV, /i/ ac = 3.5 GeV and 
Hhard = 5.0 GeV. For the error estimate we vary 
these values between 2.0 GeV < fi so ft < 3.5 GeV, 
2.0 GeV < [i Sac < 5.0 GeV and 2.5 GeV < 
[i-hard < 10.0 GeV. For the separation scale 
H se p which appears as additional parameter in 
the definition of the PS-mass we employ a value 
of fi aep = 2.0 ± 1.0 GeV. 

The analysis is performed independently in the 
pole- and PS-scheme. In figure |l| we have plotted 
the PS-mass as a function of n and the influence 
of the scale variations. The largest contribution 
to the error comes from the soft scale. Adding 
the errors from all input parameters quadrati- 
cally, our final results for the masses are 

M b = 4.984 ±0.125 GeV, 
m PS ,b (2.0 GeV) = 4.561 ± 0.112 GeV , 

m b (m b ) = 4.241 ± 0.098 GeV . (4) 

6.2. Charm mass 

As in the bottom case we use £ — 0.5. At 
this value the pole contributions still represent 
the dominant part. In principle one would like 
to choose a higher value where the theoretical 
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Figure 2. Thick solid line: central PS-mass; thin 
solid lines: mps yC for [i S oft = 1-0 and 1.25 GeV; 
dashed lines: mps, c for M/ac = 1-2 and 1.65 GeV; 
dotted lines: mps.c for Hhard = 1.4 and 2.5 GeV. 

expansions converge better. However, the con- 
tribution from the theoretical poles varies signif- 
icantly with the scales; for £ > 1 the mass de- 
pends too strongly on these variations. Thus we 
again use a range of < £ < 1. Since the per- 
turbative expansions converge more slowly than 
for the upsilon we restrict the analysis to smaller 
values of 4 < n < 7. As central values for 
our scales we have selected /i/ ac = 1.45 GeV, 
Hhard = 1-75 GeV and values of pi soft = 1-2 
GeV and fi so ft = 1-1 GeV in the pole- and PS- 
schemes respectively. For the error estimate we 
have varied the scales between 1.4 GeV < [ihard < 
2.5 GeV, 1.2 GeV < pi fac < 1.65 GeV and 
1.1(1.0) GeV < fisoft < 1.35(1.25) GeV. For the 
separation scale we choose [i sep = 1.0 ± 0.2 GeV. 
In figure || we have plotted the PS-mass and the 
corresponding error from the scales. Finally we 
obtain the masses: 

M c = 1.754 ±0.147 GeV, 
mps, c (1.0 GeV) = 1.300 ± 0.124 GeV , 

m c (m c ) = 1.188 ±0.106 GeV. (5) 

7. Conclusions 

The method of QCD sum rules is a very power- 
ful tool to extract the masses since - by the choice 
of n and £ - it can react very sensitive to thresh- 
old. Thus, large theoretical uncertainties only 
lead to a relatively small shift in the masses. We 
have tried to develop a consistent framework to 
describe the physics of the relevant energy region. 



The masses show a stable behaviour over a large 
range of n and the dominant uncertainties origi- 
nate from the threshold expansion of NRQCD. 
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